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Abstract 

A search is performed for heavy resonances decaying to two long-lived massive neu- 
tral particles, each decaying to leptons. The experimental signature is a distinctive 
topology consisting of a pair of oppositely charged leptons originating at a separated 
secondary vertex. Events were collected by the CMS detector at the LHC during 
pp collisions at a/s = 7TeV, and selected from data samples corresponding to 4.1 
(5.1) fb^^ of integrated luminosity in the electron (muon) channel. No significant ex- 
cess is observed above standard model expectations, and an upper limit is set with 
95% confidence level on the production cross section times the branching fraction to 
leptons, as a function of the long-lived massive neutral particle lifetime. 
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1 Introduction 

Several models of new physics predict the existence of massive, long-lived particles which 
could manifest themselves through their delayed decays to leptons. Such scenarios arise, for 
example, in various supersymmetric (SUSY) scenarios such as "split SUSY" [IJ or SUSY with 
very weak R-parity violation [2j, "hidden valley" models [3J, and Z' models that contain long- 
lived neutrinos H). 

This Letter presents the first search using data from the Compact Muon Solenoid (CMS) for 
massive, long-lived exotic particles X that decay to a pair of oppositely charged leptons. We 
search for events containing a pair of oppositely charged electrons or muons (dileptons) origi- 
nating from a common secondary vertex within the volume of the CMS tracker, that is signif- 
icantly transversely displaced from the event primary vertex. These leptons are assumed to 
originate from a 2-body decay of a long-lived particle, and so are required to form a narrow 
resonance in the dilepton mass spectrum. This topological signature has the potential to pro- 
vide clear evidence for physics beyond the standard model (SM). It is also very powerful in 
suppressing backgrounds from standard model processes. 

This signature is sensitive to a wide class of models. However, for the purpose of establishing a 
signal benchmark, a specific model of a long-lived, spinless, exotic particle X which has a non- 
zero branching fraction to dileptons is used. In this particular model, the X is pair-produced 
in the decay of a (non-SM) Higgs boson, i.e. H*^ — > 2X, X — > [jS], where the Higgs boson 
is produced through gluon-gluon fusion. This model predicts up to two displaced dilepton 
vertices in the tracking volume per event. 

The DO Collaboration has performed searches for leptons from delayed decays in its tracker 
volume in [3, but these searches are sensitive to a much smaller kinematic phase space region 
than CMS. The ATLAS Collaboration has performed searches that are sensitive to decay lengths 
up to about 20 m by exploiting the ATLAS muon spectrometer [81 19|, using different decay 
channels from those considered in this Letter. 

2 The CMS detector 

The central feature of the CMS apparatus fTO\ is a superconducting solenoid of 6 m internal 
diameter providing an axial field of 3.8 T. Within the field volume are the silicon pixel and strip 
tracker, the lead-tungstate crystal electromagnetic calorimeter (ECAL), and the brass/ scintillator 
hadron calorimeter. Muons are identified in gas-ionisation detectors embedded in the steel 
magnetic-flux return yoke of the solenoid. 

The silicon tracker is composed of pixel detectors (three barrel layers and two forward disks 
on either end of the detector) surrounded by strip detectors (ten barrel layers plus three in- 
ner disks and nine forward disks at each end of the detector). The tracker covers the pseu- 
dorapidity range \t]\ < 2.5, where t] = — ln[tan(0/2)] and 6 is the polar angle with respect 
to the anticlockwise-beam direction. All tracker layers provide two-dimensional hit position 
measurements, but only the pixel tracker and a subset of the strip tracker layers provide three- 
dimensional hit position measurements. Owing to the strong magnetic field and the high gran- 
ularity of the silicon tracker, promptly produced charged particles with transverse momentum 
Pj = lOOGeV/c are reconstructed with a resolution in pj of ?»1.5% and in transverse impact 
parameter do of ~15 f/m. The track reconstruction algorithms are able to reconstruct displaced 
tracks with transverse impact parameters up to «25 cm from particles decaying up to ~50 cm 
from the beam line. The performance of the track reconstruction algorithms has been studied 
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with data (TTl .The silicon tracker is also used to reconstruct the primary vertex position with a 
precision of cr^j ^ 20 }im in each dimension. 

The ECAL consists of nearly 76 000 lead tungstate crystals, which provide coverage in pseudo- 
rapidity \f]\ < 3. Muons are measured in the pseudorapidity range \rj\ < 2.4 with detection 
planes based on one of three technologies: drift tubes in the barrel region, cathode strip cham- 
bers in the endcaps, and resistive plate chambers in the barrel and endcaps. 

The first level of the CMS trigger system, composed of custom hardware processors, selects 
events of interest using information from the calorimeters and the muon detectors. A high- 
level trigger processor farm then employs the full event information to further decrease the 
event rate. 

3 Data and Monte Carlo simulation samples 

For this analysis, pp collision data at a centre-of-mass energy of 7TeV corresponding to an 
integrated luminosity of 4.1 ± 0.1 (5.1 ± 0.1) fb^^ are used in electron (muon) channels. (The 
lower electron luminosity is due to the fact that not all data was taken with triggers suitable for 
this analysis.) 

For the electron channel, these data are collected with a trigger that requires two clustered en- 
ergy deposits in the ECAL, each with transverse energy Ej > 38 GeV. For the muon channel, 
the trigger requires two muons, each reconstructed without using any primary vertex con- 
straint and having pj > 30 GeV/c. The tracker information is not used in either trigger. 

Signal samples are generated using PYTHIA V6.424 [12 J to simulate H'^ production through 
gluon fusion (gg — > H*^). Subsequently the H" is forced to decay to two long-lived spin 
exotic particles (H*^ XX), which then decay to dileptons (X i^(^), where £ represents 
either a muon or an electron. Several samples with different combinations of masses (M^o 
= 125, 200, 400, 1000 GeV/c^) and X boson masses (Mx = 20, 50, 150, 350 GeV/c^) are generated. 
The lifetimes of X bosons used in these samples are chosen to give a mean transverse decay 
length of approximately 20 cm in the laboratory frame. Several simulated background samples 
generated with pythia are used, corresponding to tt, Z/ 7 — > £^i^ (including jets), W/Z boson 
pair production with leptonic decays, and QCD multijet events. The contribution from single 
W + jet production is negligible. In all the samples, the response of the detector is simulated in 
detail using Geant4 ||T3| . The samples are then processed through the trigger emulation and 
event reconstruction chain of the CMS experiment. 

4 Event reconstruction and selection 

Events are required to contain a primary vertex with at least four associated tracks whose po- 
sition is displaced from the nominal interaction point by no more than 2 cm in the direction 
transverse to the beam and no more than 24 cm in the direction along the beam. Furthermore, 
to reject events produced by the interaction of beam-related protons with the LHC collimators, 
the fraction of tracks classified as "high purity", as defined in Ref. |ll| , must exceed 25% in any 
event with at least 10 tracks. 

The selection of lepton candidates from displaced secondary vertices begins by searching for 
high-purity tracks with transverse momenta pj > 41 GeV/c (33 GeV/c) for the electron (muon) 
channel. These criteria are slightly higher than the corresponding trigger thresholds, to min- 
imise dependence on the trigger inefficiency in the pj turn-on region. The tracks are required 



4.1 Selection efficiency 



3 



to have pseudorapidity < 2, as the efficiency for finding tracks from displaced secondary 
vertices falls off at large To reject promptly produced particles, the tracks must have a 
transverse impact parameter significance with respect to the beam line of |do/c^d| > 3 (2) in the 
electron (muon) channel. Because bremsstrahlung in the material of the tracker significantly 
affects the reconstruction of electrons, the electron channel selection criteria are more restric- 
tive. Tracks are considered to be identified as leptons if AR = ^ (A?/)^ + (A(/))^ is less than 
0.1. Here, A?/ and A^ are the differences between the track and a lepton trigger object in pseu- 
dorapidity and ^, the azimuthal angle about the anticlockwise-beam direction. Standard CMS 
offline lepton identification algorithms are not applied, since they are inefficient for leptons 
from highly displaced vertices. However these algorithms are not needed to suppress the very 
low backgrounds present in this analysis. For the electron channel, the energy of the electron 
is estimated from a deposit in the ECAL that is near the reconstructed track trajectory. 

The X boson candidates are formed from pairs of oppositely-charged lepton candidates. The 
two corresponding tracks are fitted to a common vertex, which is required to have a chisquare 
per degree of freedom ;^-^/dof < 5. For events in the electron (muon) channel, the vertex 
must lie at a distance of more than 8 (5) standard deviations from the primary vertex in the 
transverse plane. If either track has more than one hit closer to the centre of CMS than their 
common vertex, the event is rejected. 

Both lepton candidates are required to be isolated, to reject background from jets. A hollow 
isolation cone of radius 0.03 < Ai? < 0.3 is constructed around each candidate. Within this 
isolation cone, the Ya Pt of tracks with pj > 1 GeV/ c, excluding the other lepton candidate, 
must be less than 4 GeV/ c. This requirement has very little effect on the signal efficiency, which 
is relatively insensitive to the number of primary vertices in each event. According to simu- 
lation, the mean J^pj in the isolation cone increases from 0.6 to 1.2GeV/c as the niunber of 
additional primary vertices per event increases from to 20. 

Cosmic ray muons may be reconstructed as back-to-back tracks. To reject them, a requirement 
of cos (a) > —0.95 is applied, where a is the opening angle between the two tracks. Background 
from misidentified leptons is reduced by requiring that the two lepton candidates are not both 
matched to the same trigger object. Additional background rejection is achieved by requiring 
that, projected into the plane perpendicular to the beam line, the reconstructed momentum 
vector of the X boson candidate is collinear with the vector from the primary vertex to the 
secondary vertex. The coUinearity angle is required to be less than 0.8 (0.2) radians in the 
electron (muon) channel. Owing to the difficulty of modelling the trigger efficiency for closely 
spaced muon pairs, the two tracks in muon channel candidates must be separated by AR > 0.2. 
To eliminate background from ]/xp and Y decays and from 7 conversions, X boson candidates 
are required to have dilepton invariant masses larger than 15 GeV/c^. If more than one X boson 
candidate is identified in a given event, all the selected candidates are retained. 

4.1 Selection efficiency 

The efficiency of the full set of selection criteria for X £^£^ decays is defined on a per- 
candidate basis, not on a per-event basis. It is determined from the simulation, separately 
for the electron and muon channels, for two different classes of events: first for H° — > XX 
events in which only one long-lived exotic particle decays to the chosen lepton species, defining 
efficiency ei, and second for events in which both long-lived exotic particles decay to chosen 
lepton pairs, defining efficiency £2- The efficiencies ei and £2 are usually almost identical, 
indicating that the efficiency to select an X boson candidate is not strongly affected by whether 
or not the second X boson in the event decays to the same lepton channel. The only exception is 
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for the muon channel in the case of small Mx/ M^o, where the dimuon trigger is inefficient for 
the two nearly collinear muons from the decay of the same X boson, but the trigger requirement 
can still be satisfied by muons from separate X bosons. The efficiencies are estimated for a 
range of X boson lifetimes, corresponding to mean transverse decay lengths of ~ 0.7 - 600 cm, 
by reweighting the simulated signal events. The maximum efficiency (for M^p = lOOOGeV/c^, 
Mx = 150 GeV/c^, ct = 1 cm) is approximately 34% (52%) in the electron (muon) channel, but 
becomes significantly smaller at lower H*^ masses or longer lifetimes. 

5 Background estimation and modelling 

An interpretation of the observed dilepton mass spectrum requires an estimate of the back- 
ground normalisation and a parametrisation of the background shape as a function of Mx- 

5.1 Background normalisation 

The number of background events passing all the selection criteria for X boson candidates, is 
estimated from simulated samples using the distribution of the transverse decay length signif- 
icance Lxi// Cxy We parameterize this distribution with the sum of two falling exponentials. By 
integrating the fitted curve over the signal region, defined by hxyl o'xy > 8 (5) for the electron 
(muon) channel, an estimate of the mean total background in the mass spectrum is obtained. 
The estimate gives 1.4^J 2 (0.02^q q2) candidates in the electron (muon) channel. This estimate 
of the mean total background is used to derive the results in Section|7| To verify that the simula- 
tion correctly describes the Lxy/ (Txy spectrum, data and simulation distributions are compared 
in Fig. [l| after relaxing the lifetime-related selection requirements to increase the number of 
events and ensure that the plots are dominated by background. Specifically, the threshold on 
the decay length significance Lxy/ (Txy was decreased to 3, for both the electron and muon chan- 
nels, and the requirement on the individual lepton do/ cT'd was removed entirely. Figure [l] also 
shows that the main background to this search consists of prompt dileptons that have been 
reconstructed with large decay length significance. 
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Figure 1: The transverse decay length significance of the candidates for the dielectron (left) and 
dimuon (right) channels with loosened cuts in data and simulation. The vertical dashed line 
indicates the selection requirement used for signal events. There are no simulated QCD or tt 
events passing these selection requirements, so they are omitted. 
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5.2 Background shape 

An estimate of the background shape can be obtained from the Mx distribution of a back- 
ground sample. However, after applying all selection requirements there are too few events 
to measure its shape accurately. Since the dilepton mass distribution and lifetime-related vari- 
ables are only weakly correlated in simulated background candidates, the shape of the mass 
distribution is instead obtained by fitting a parameterized function to data samples with the 
lifetime-related selection requirements removed. Namely, no selections are made on the indi- 
vidual lepton do/ dd, the transverse decay length significance Lxy/ dxy, or the collinearity angle 
A(p. Figure |2] shows the results of these fits to the electron and muon data samples. The back- 
ground is represented by the sum of two functions: a Breit-Wigner function, to represent the 
Z resonance, multiplied by a Gaussian error function, to approximate the effect on the selec- 
tion efficiency of the lepton pj thresholds; and an additional function, given by an exponential 
term, to represent the non-Z background. The fits give the fraction of the background from the 
Z as 0.985 ± 0.002 (0.994 ± 0.001) for the electron (muon) channels. 
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Figure 2: Distribution of the dilepton mass and the fitted shape in a data sample with lifetime- 
related selection requirements removed, shown for the electron (left) and muon (right) chan- 
nels. The shape used is that of a Breit-Wigner distribution times a turn-on function, added to 
an exponential term. 



6 Systematic uncertainties 

The primary systematic uncertainty comes from the efficiency in detecting and reconstructing 
signal events. This uncertainty derives from uncertainties in the efficiency of reconstructing 
tracks from displaced vertices, the trigger efficiency, the modelling of pileup in the simulation, 
the parton distribution function sets, the renormalisation and factorisation scales used in gener- 
ating simulated events, and the effect of higher order QCD corrections. In addition, systematic 
uncertainties in the integrated luminosity, and the background estimate are considered. 

Table [l] summarises the sources of systematic uncertainty affecting the signal efficiency. The 
relative uncertainty in the luminosity is taken to be 2.2% [14J. 

Varying the modelling of the pileup within its estimated uncertainties yields a relative change 
in the signal selection efficiency of less than 2%, irrespective of the mass point chosen. The 
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Table 1: Systematic uncertainties affecting the signal efficiency over the range of M^o and Mx 
values considered. In all cases, the uncertainty specified is a relative uncertainty. Note that the 
NLO uncertainty is only evaluated for the Mj^o = 125 GeV/c^ case. The relative uncertainty in 
the luminosity is taken to be 2.2%. 



Source 


Uncertainty 


Pileup modelling 


2% 


Parton distribution functions 


<1% 


Renormalisation and factorisation scales 


<0.5% 


Tracking efficiency 


20% 


Trigger efficiency 


2.6% (e), 11% {^) 


NLO effects (Mho = 125 GeV/c^ only) 


4-12% 



relative uncertainty due to parton distribution functions (PDF) is studied using the PDF4LHC 
procedure [15| and is less than 1% for all mass points. The dependence of the acceptance on 
the choice of the renormalisation and factorisation scales, which are chosen to be equal and are 
varied by factors of 0.5 and 2, is found to be well below 0.5%. These uncertainties are applied 
in the cross-section limit calculation. 

6.1 Track-finding efficiency 

Understanding the efficiency to find a track as a function of its impact parameter is a crucial 
aspect of the analysis. Two methods are used to assess if the efficiency to reconstruct displaced 
tracks is correctly modelled by the simulation. 

The first method consists of a direct measurement of the efficiency to reconstruct isolated tracks, 
using cosmic ray muons. Events are selected from dedicated runs with no beam activity and the 
cosmic ray muons are reconstructed by combining the hits in the muon chambers from opposite 
halves of CMS. The efficiency to reconstruct a tracker track associated with a cosmic ray muon, 
as a function of the transverse and longitudinal impact parameters, is shown in Fig.|3] We focus 
principally on the region |£fo| < 20 cm, since in simulated signal, the reconstructed tracks from 
displaced vertices lie predominantly in this region. Data and simulation agree within 10% in 
this region, so the corresponding relative systematic uncertainty in the efficiency to reconstruct 
dilepton candidates is thus taken to be 20%, as there are two tracks per candidate. This method 
does not explicitly measure tracking efficiency for dielectron candidates. However, simulation 
studies indicate that the electron efficiency is only about 10% smaller than the muon efficiency, 
where the difference can be attributed to bremsstrahlung in the electron case. The material 
budget in the tracker is modelled in simulation with an accuracy better than 10% ||T6| . It is 
assumed that the difference in tracking efficiency between electrons and muons is modelled 
with similar precision. The difference can therefore be neglected in comparison with the much 
larger systematic uncertainty in the efficiency to reconstruct dilepton candidates. 

To determine if the simulation properly describes the tracking efficiency in the presence of a 
high density of hits, a second method is used. Single charged particles are simulated at various 
production points throughout the tracker volume. These particles are then embedded in both 
data and simulated data events, and the difference in efficiency to reconstruct the particles is 
compared. The relative difference between data and simulation is less than 1%. 

6.2 Trigger efficiency measurement 

The trigger efficiency is measured using the "tag-and-probe" method flT]. The decays of 
Z bosons to dileptons are reconstructed in data collected with single-lepton triggers. They 
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Figure 3: Efficiency of the tracker to find a track, given a cosmic ray muon reconstructed in the 
muon chambers, as a function of the transverse (left) and longitudinal (right) impact parame- 
ters (with respect to the nominal interaction point of CMS). The efficiency is plotted in bins of 
2 cm width. For the left plot, the longitudinal impact parameter |zo| is required to be less than 
10 cm, and for the right plot, the transverse impact parameter |do | must be less than 4 cm. 

are used to measure the efficiency for a lepton to pass one leg of the dilepton triggers used in 
this analysis. The dilepton trigger efficiency is then obtained as the square of this single-lepton 
efficiency, which assumes that there is no correlation in efficiency between the two leptons. 
This is generally a good assumption except in the muon channel, where dimuons separated by 
AR < 0.2 must be excluded because the trigger is inefficient for closely spaced dimuons. The 
systematic uncertainty associated with the trigger efficiency is evaluated by taking the differ- 
ence between the estimates of the efficiency from data and simulation, yielding a total relative 
uncertainty of 2.6% for the electron channel and 11% for the muon channel. 

6.3 Effect of higher-order QCD corrections 

For Myf = 125GeV/c^, the leptons from the X boson decay have a combined efficiency of 
only a few percent for passing the lepton pj requirements. For this reason the signal efficiency 
at this mass is sensitive to the modelling of the Higgs px spectrum, which may in turn be 
influenced by higher order QCD corrections. To study this effect, we reweight the LO H'^ pj 
spectrum from our signal sample to match the corresponding Higgs pj spectrum evaluated at 
NLO. For Mho = 125GeV/c2 and Mx = 20 (50) GeV/c^ the signal efficiency changes by 4% 
(12%). This change is taken as an additional systematic uncertainty in the efficiency for the 
M^o = 125GeV/c^ case. For larger H*^ masses, the corresponding systematic uncertainty is 
below 0.5%, and hence neglected. 

6.4 Bacltground uncertainty 

The systematic uncertainty for the background normalisation is taken from the uncertainty on 
the background fits described in Section 5.1 where these fits yield 1.4+^ 2 (0.02^0 02) candidates 
in the electron (muon) channel. In addition, a cross-check is made by repeating the fit described 
in Section 5.1 but restricting the fit to the background region Lxy/ (Txy < 8 (5) in the electron 



(muon) channel and then extrapolating it into the signal region. This cross-check is performed 
for both simulated events and data. Additionally, the background normalisation is estimated 
by simply counting the number of candidates in the simulated background passing all selection 
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criteria. All methods give consistent results. As explained in Section the uncertainty in the 
background normalisation actually affects the observed limits only for X boson masses close to 
the Z resonance. 

To determine the sensitivity of the limits to variations in the assumed background shape, de- 



scribed in Section 5.2 other functional forms are also tried. The resulting effect on the limits, 
which is taken as a systematic uncertainty, is negligibly small. The fit for the background shape 
is repeated using simulated events and consistent results are obtained. As an additional check 
that the shape is not strongly influenced by the lifetime-related requirements, the mass distri- 
bution obtained by relaxing the lifetime-related selection requirements is also fitted, in both 
data and simulated events. In all cases, the resulting change in the limits is negligible. 



7 Results 

After all selection requirements are applied, no candidates survive in the muon channel, consis- 
tent with the expected mean number of 0.02^q I^i the electron channel, a total of 4 candidates 
remain, which is also consistent with the expected mean number of 1.4^| However, 2 of these 
candidates fall in the region above the Z mass, where the estimated background is lower. The 
observed candidates all have transverse decay lengths less than 0.3 cm. 

We set 95% confidence level (CL) upper limits on the signal process using the CLg method 
||T8l[l9J , which makes use of an unbinned likelihood fit to the dilepton mass spectrum. 

This fit to the mass spectrum uses the following functions: 

• A Gaussian signal function to represent the signal's mass distribution. For each 
H*^ mass, the mass resolution used in the Gaussian is obtained from the simulated 
signal samples, as a function of the X boson mass and lifetime, interpolating between 
the generated X boson masses with a smooth curve when necessary. 

• The sum of two background functions, one distribution representing the background 
from the Z peak and another more slowly varying distribution representing the non- 
Z background. These functions are obtained as described in Section |5] 

The limit calculation takes into account the systematic uncertainties described in Section |6] by 
introducing a nuisance parameter for each uncertainty, marginalized by a log-normal prior 
distribution. In particular, the normalisation of the Z (non-Z) background can be constrained by 
the a priori estimate of the total background normalisation presented in Section [SH^ multiplied 
by the estimated Z (non-Z) background fraction from Section [5^ However, the fits to the mass 
spectrum performed as part of the limit calculation strongly constrain the normalisation of the 
non-Z background, even though the normalisation of a signal is unknown. As a result, it is not 
necessary to use the a priori estimate of the non-Z background normalisation when calculating 
the observed limits. In the case of the Z background, the a priori estimate of the normalisation 
is used, but it affects only the limits for X bosons whose mass is close to that of the Z. To 
calculate expected limits, one must have a prediction of the background normalisation, which 
is taken from Section |5l 

As a first step, upper limits are placed on the mean number Nx of X bosons that could pass 
the selection requirements, as a function of the X boson mass. The resulting upper limits on 
Nx at 95% CL for the electron and muon channels are presented in Fig. |4j These limits are 
independent of the particular model assumed for X boson production, except for the mass 
resolution assumed for the signal, which affects the width of the resulting peaks in the observed 
limit. The mass resolution used for these limits is derived from the Monte Carlo simulation for 
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the hypothesis Mh = lOOOGeV/c^, which has the largest mass resolution of the signal points 
studied and hence yields the most conservative limits. The limits on Nx are close to 3.0 at 
most masses, as one would expect from Poisson statistics with zero observed signal, but are 
larger at mass points near the masses of dilepton candidates observed in the data. Since the 
fitted background levels under the signal peak are extremely small, the limits are not expected 
to depend on the background shape, and indeed using the alternatives described in Section |5] 
give negligible changes in the results. This figure also shows the 95% CL expected limit band. 
Except near the Z resonance, the a priori predictions of the background normalisation are very 
small, so the expected limit is close to 3.0 and the expected limit band is extremely narrow; the 
median value of the expected limit is in fact equal to 3.0 ever3rwhere. 
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Figure 4: The 95% CL upper limits on the mean number of X bosons that could pass the selec- 
tion requirements in the electron (muon) channels are shown in the left (right) plot. A yellow 
shaded band shows the 95% quantile for the expected limits, but is almost entirely hidden by 
the observed limit curves. 

The expected number of signal dilepton candidates passing the selection cuts can be expressed 
as: 



Nx = 2Cei(rB[l + B{e2/ei-l)] 



4.1 



(1) 



cr is the 



where C is the integrated luminosity, e(i,2) the efficiencies defined in Section 
production cross section of the heavy resonance decaying to X X, and B is the branching fraction 
for the decay X — ?> i^£^ . This expression takes into account that either one or both X bosons 



in an event may decay to the chosen lepton species, and that, as mentioned in Section 4.1 the 
efficiency to select such an X boson is slightly different in the two cases. Using this equation, 
the likelihood function can be expressed in terms of crB, thus allowing upper limits to be placed 
on this quantity. Since Nx in Eq. ([T]) depends not only on uB, but also on B, the upper limits 
depend on the assumed value of B. However, the factor ( 62/61 — 1) is in practice always 
positive or very small. Hence if one assumes infinitesimally small B when calculating the 
limits on aB, such that the factor in square brackets in Eq. ^ is equal to 1, the resulting limits 
will be valid, and in some cases conservative, for any value of B. 

For each combination of the H^ and X boson masses that are modelled, and for a range of X 
boson lifetimes, the 95% CL upper limits on aB are calculated. The observed limits are shown 
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in Figs, jsj-jzl (No results are shown for M^o < 200 GeV/ in the electron channel, since the high 
trigger thresholds result in a very low signal efficiency.) Note that for the muon channel in the 
Mho = 1000GeV/c2, Mx = 20GeV/c2 case, the efficiency is significantly reduced because the 
muons are produced very close together, which causes trigger inefficiencies. Since the observed 
dilepton candidates do not have masses close to those of the X bosons considered in these plots, 
they have no effect on the limits. The bands show the 95% quantile for the expected limits. 
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Figure 5: The 95% CL upper limits on aB for the electron (left) and muon channel (right) for a 
H*^ mass of 1000 GeV/c^. Narrow yellow shaded bands show the 95% quantiles for the expected 
limits. 

For H*^ or X boson masses other than those plotted in Figs, [s] - 17| exact limits are not com- 
puted, since no simulated signal samples are available with which to determine the signal 
selection efficiency. However, since the observed limits appear to be monotonic functions of 
the H*^ and X boson masses, one can infer approximate limits for other masses, provided the 
latter lie within the range of those shown in the figures. For example, for Mfjo = lOOOGeV/c^, 
it should be safe to assume that the limits for Mx = 170GeV/c^ would be at least as good as 
the weaker of the limits for Mx = 150GeV/c^ and Mx = 350GeV/c^. However, for X bosons 
that are close in mass to the candidates seen in data, the limits would be worse than this. For 
these particular masses. Fig. |4] gives an indication of the factors by which the limits would be 
degraded. 

In an alternative signal model, in which H*^ — ?► XY, where the Y boson does not decay to dilep- 
tons, the expected number of signal candidates Nx passing the selection cuts can be expressed 
as Nx = Ce[crB, where e[ is the efficiency to select signal candidates, in this scenario. If the X 
and Y bosons have identical masses, then = ei, and comparison with Eq. (jljl shows that the 
limits would be a factor of 2 worse than those presented in Figs.|5]-[7[ 

The limits quoted above are for H*^ bosons produced through gluon-gluon fusion. If the H*^ 
bosons were instead produced by the sum of all standard model production mechanisms, their 
momentum spectra would be different. For M^o = 125GeV/c^, the selection efficiency would 
then be larger by a factor of approximately 1.18 (1.08) for Mx = 20 (50) GeV/c^ and there would 
be a corresponding improvement in the limits. 

If the initial resonance were a Z' with spin 1 instead of a H*^ boson, the acceptance for a dilepton 
pair to pass the pj and rapidity selection cuts would be slightly different, mainly because Z' 
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mass of 400 GeV/c^. Narrow yellow shaded bands show the 95% quantiles for the expected 
limits. 
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are produced by qq annihilation, whilst H are produced predominantly through gluon-gluon 
fusion. If the Z' is to decay to a pair of long-lived spin particles, these cannot be identical if 
CP is to be conserved, but they are assumed in what follows to have equal mass. Their angular 
distributions would differ depending on whether they were fundamental spin bosons or 
spin 2 'hidden valley' quarks that have hadronised in the dark sector into spin bosons ||3|. 
Studies with simulated events show that the change in the acceptance for these two Z' models 
relative to the original H*^ model is less than 3% for a Z' mass of 1000 GeV/c^, less than 11% for 
a 400 GeV/ Z' mass, and below 25% (in the muon channel) for a 200 GeV/ Z' mass. Thus the 
limits quoted above on the H*^ model will be approximately valid for these two Z' models as 
well, for Z' masses of 200 GeV/c^ or more. 

8 Summary 

A search for long-lived neutral particles, X, produced in pp collisions at ^/s = 7 TeV and decay- 
ing to either e+e^ or has been performed. In the e+e^ channel 4 candidates are observed, 
of which 2 are in the Z mass region, and in the ji^ji^ channel no candidates are observed. These 
results are consistent with standard model expectations and are used to derive upper limits on 
the product of cross section times branching fraction for a Higgs boson, in the mass range 200 - 
1000 GeV/c-^, decaying into a pair of X bosons, in the mass range 20-350 GeV, which each decay 
to e+e^ and ■ The limits are typically in the range 0.7-10 fb, for X bosons with lifetimes in 
the range 0.1 < ct < 200 cm. For a Higgs mass of 125GeV/c^, the corresponding limits are in 
the range 10-100 fb. These are the most stringent limits in these channels to date. 
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